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Geometry induced entanglement transitions in nanostructures 

J. p. Coe,i'E3S. Abdullah, i'E3 and I. D'Amicoi'El 

Department of Physics, University of York, York YOlO 5DD, UK 

We model quantum dot nanostructures using a one-dimensional system of two interacting electrons. We show 
that strong and rapid variations may be induced in the spatial entanglement by varying the nanostructure 
geometry. We investigate the position-space information entropy as an indicator of the entanglement in this 
system. We also consider the expectation value of the Coulomb interaction and the ratio of this expectation to 
the expectation of the confining potential and their link to the entanglement. We look at the first derivative of 
the entanglement and the position-space information entropy to infer information about a possible quantum 
phase transition. 
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I. INTRODUCTION 

For a long time entanglement had been thought of as 
a failure of quantum mechanics or at best a quantum cu- 
riosity, but it is now considered a physical resource which 
could feasibly allow quantum information processors to 
overcome some of the limitations of their classical ana- 
logues. 

Semiconductor quantum dots have been put forward 
as possible components of solid state quantum informa- 
tion processing devices^"— . One of the possible advan- 
tages of the use of quantum dots is that the param- 
eters of the system may be tuned, allowing the prop- 
erties of semiconductor nanostructures to be tailored^. 
The seemingly inexorable progress of technology ap- 
pears to promise advanced engineering of quantum dot- 
based structures, thus leading to the fabrication of cou- 
pled and scalable quantum dot systems. To use quan- 
tum dot devices for quantum computing necessitates the 
ability to generate and manipulate entanglement within 
these structures. In quantum dot systems this entan- 
glement could be controlled through externally applied 
electro-magnetic fields^ or by varying the parameters of 
the nanostructure^. In this contribution, we show how 
changing the geometry of the confinement potential of 
single, core-shell and double quantum dot structures can 
influence the spatial entanglement^ between two elec- 
trons trapped within the nanostructure; such a property 
could possibly be exploited to design nanostructures ac- 
cording to the level of entanglement needed for a specific 
application. We also investigate the use of the position- 
space information entropy as an indicator of the entan- 
glement for this system. Our results show a potential 
quantum phase transition induced by varying the nanos- 
tructure confinement geometry. A quantum phase tran- 
sition is characterized by a nonanalyticity in the ground 
state energy and could be detected by the nonanalyticity 
of the entanglement measure^ which we seem to observe. 



The transition we observe could be created experimen- 
tally by changing the polarity of a gate positioned over 
the core-to-barrier region within a gate-defined quantum 
dot. 



II. THE MODEL 

We model the quantum dot system with a one dimen- 
sional Hamiltonian describing two interacting electrons 
constrained by a two-center power-exponential potential 



2dx^ 



U{xi,X2) 



(1) 



where we employ (effective) atomic units. Here we model 
the Coulomb repulsion via the contact interaction^""— 
U{xi,X2) — S{xi — X2) while the two-center power- 
exponential potential^^ V{x) is given by 

V{x) = -Vo {exp[-(|a; + d\/Rf] + exp[-(|a; - d\/R)P]} . 

(2) 

This potential allows us to access a range of nanostruc- 
ture geometries: single dots, double dots and core-shell 
dots (well-within-a-well)§-. We may move continuously 
from one structure to the next by varying the potential 
range R. For fixed potential well depth Vq and well cen- 
ter distance 2d, the parameter p allows us to control the 
'hardness' of the potential. This, for example, allows 
us to consider quantum dots with little or no intermix- 
ing between the different materials using a large p value 
(hard potential, p=200) or, by using a smaller p value, 
p < 7, smoother ('softer') potentials corresponding to the 
experimental situation where there is substantial mixing 
at the interface between materials. We set p = 200 so as 
to have almost square wells and thereby allow strong and 
very fast transitions in the entanglement to be observed. 
We solve the system using 'exact' diagonalization based 
on a single-particle harmonic oscillator basis set^. 



III. ENTANGLEMENT 
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As the spin-part of the ground state wavefunction is 
constant (in fact a singlet: a maximally entangled qubit 



2 



state) we investigate the spatial entanglement^. We mea- 
sure the entanglement using the linear entropy of the re- 
duced density matrix 

L = Tr{p,^d - p?ed) = 1 - Trp^.^^. 

Here the reduced density matrix is prcd — Trg \^) {ipl 
and the subsystem B corresponds to one of the electrons. 
It has been shown that L may be used as a valid measure 
of entanglement when quantifying the particle-particle 
entanglement!^^ Importantly L is much easier to cal- 
culate than S, the von Neumann entropy of the reduced 
density matri»i^, when a system has a very large number 
of degrees of freedom. This is the case for the spatial en- 
tanglement where the degrees of freedom are continuous. 
In this case the reduced density matrix is given by 

with the many-body wave-function, 

p'^^^{xi,X2) = J Picdixi,X3)p^cdix3,X2)dX3, 

and 



For a maximally entangled system with continuous de- 
grees of freedom L would be equal to unity. 

We also compare L with the position-space information 
entropy of the density 



Sn 



n{x) log n{x)dx, 



(3) 



where n{x) is the particle density. This has been used 
to study the entanglement of the Moshinsky atom^** and 
may also be thought of as a zeroth-order approximation 
to S, where the off-diagonal terms of the reduced density 
matrix are set to zeroi^ 



IV. RESULTS 

Our results show that the spatial entanglement can in- 
deed be successfully manipulated through the geometry 
of the confining potential: in Fig. [TJ upper panel, we 
show that, by manipulating the geometry of the nanos- 
tructure, we can create entanglement, as quantified by 
L, from almost zero to 0.5 (the maximum for a system of 
two qubits). For R ^ SO the potential is a well-within- 
a-well for which the inner well becomes narrower with 
decreasing R (see inset of Fig. [T]) and the entanglement 
decreases until it is almost zero. In this case the elec- 
trons are now essentially described by a product state, 
the wave function dominated by the very strongly con- 
fining potential. This inner well then inverts to become 
a barrier and gives a double well structure which has an 
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FIG. 1. Upper panel: The linear entropy of the reduced den- 
sity matrix L and the position-space information entropy of 
the density S„ scaled by 0.1. Lower Panel: The partial deriva- 
tive with respect to i? of L and Sn. Inset: potential shape for 
different values of R, as labeled. 



entanglement of 0.5. This entanglement is relatively 
robust: as R decreases further, the distance and width 
of the double wells change (see inset of Fig. [T]), but the 
entanglement does not. The sharp transition from al- 
most zero entanglement to the maximal entanglement 
for a qubit system occurring at i? ~ 8 in Fig. [T] means 
that we have moved from an essentially non-interacting 
non-entangled product state - due to the very strongly 
confining inner potential well in the core-shell QD - to 
a highly entangled state where the wave-function splits 
between the two wells. We find that a discontinuity in 
the derivative (lower panel of Fig. [1]) of the entanglement 
measure with respect to the potential range R seems to 
begin to appear at this point. This discontinuity suggests 
a possible quantum phase transitionii 

When comparing Sn, the position-space information 
entropy (Eq. [3]), to the linear entropy L, we find that, 
for R greater than the transition point, Sn has a behav- 
ior qualitatively similar to L. However for lower R it 
does not display a constant value but instead sharply de- 
creases. This confirms that Sn cannot be considered as 
a measure, nor an indicator of entanglement. It is in- 
teresting to notice though that the sharp entanglement 
transition at i? « 8 is captured by the position-space 
information entropy Sn- 

Our results suggest that the entanglement has an in- 
verse correlation with the Coulomb interaction energy 
for this system, shown in the upper panel of Fig. [2] The 
minimum of the interaction energy corresponds in fact to 
the maximum of the entanglement. This seems counter- 
intuitive as a non-zero spatial entanglement requires an 
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FIG. 2. Upper panel: The expectation value of the Coulomb 
interaction. Lower panel: The ratio of expectations of the 
Coulomb interaction and potential energy where the potential 
is shifted by its lowest value Vq. 

interaction betvifccn the particles; however we note that 
while the Coulomb interaction strength is fixed, its ex- 
pectation value will be very low when the probability of 
finding the particles in the same region is almost zero. At 
the transition (i? « 8) the interaction energy falls from 
its maximum value - where both electrons are confined in 
the inner well of a core-shell dot and there is almost zero 
entanglement- to almost zero interaction energy, when 
the electrons are confined in a double well system and 
their repulsion means they are unlikely to be found in 
the same well. 

In previous work^ we suggested that for the three- 
dimensional Hooke's atom the ratio of the Coulomb in- 
teraction energy to the potential energy could be used to 
indicate the entanglement. This was for a harmonic po- 
tential which always had zero as its lowest value. Hence 
to allow a fair comparison we shift the potential of the 
system under investigation by its lowest value Vq, and 
plot < U > / < - Vb > in the lower panel of Fig. H 
Due to the finite numerical precision in the difference 
< V > — Vb, we have to limit the data to i? < 20. We see 
that, although < U > / < V — Vq > indicates the trend of 
the entanglement for i? > 8 it does not pick up the large 
increase in entanglement for smaller R. We suggest this 



could be a consequence of the system topology and merits 
further investigation. In addition this quantity does not 
experience a sharp change as the entanglement rapidly 
shifts from its minimum to its maximum, rather it ap- 
pears to still be a smooth function in the range R ^ 8. 

V. CONCLUSIONS 

The results obtained for our model system results sug- 
gest that strong and rapid variations in the spatial en- 
tanglement in nanostructures may be induced by vary- 
ing the potential geometry. We found that the maxi- 
mum entanglement we could reach was an entanglement 
of L = 0.5 (the entanglement of two maximally-entangled 
qubits). This occurred when the nanostructure was com- 
posed of a double well and this result was robust against 
the shape and distance of the two wells. We have seen 
that the Coulomb interaction energy appears to indicate 
the entanglement through an inverse relationship; how- 
ever neither the ratio of the interaction energy to po- 
tential energy nor the position-space information entropy 
could be considered as an entanglement indicator for this 
one-dimensional system. Our numerical results suggested 
that there might be a discontinuity in the first derivative 
of the entanglement with respect to R for even harder 
confining potential {p — ?> oo). This was suggestive of a 
potential quantum phase transition and corresponded to 
the transition between minimum and maximum entan- 
glement. Interestingly, a similar feature appeared also 
in the first derivative of the position-space information 
entropy with respect to R. 
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